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1. Introduction 

The Dyson-Schwinger equation (DSE) of a gauge boson was recently revisited [1] in the con¬ 
text of the Kugo-Ojima (KO) confinement criterion [2]. It was found that the term which saturates 
the rhs. of the DSE for 0 depends on the phase of a gauge theory [1]. In the Higgs phase, for 
example, physical states saturate the transverse DSE, while in the confining phase only unphysical 
degrees of freedom confribufe fo fhe safurafion of fhe DSE. Inferesfingly, fhis corollary fo fhe KO 
criferion is nol only frue for linear covarianf gauges, buf can also be applied fo other gauges and 
models like maximal Abelian gauge, non-covariant Coulomb gauge and in the Gribov-Zwanziger 
theory [1]. 

In an earlier attempt [3] we have tried to verify the KO criterion directly on the lattice. There 
however we could not confirm fhe desired resulf for fhe KO funcfion: u{p) —)> — 1 for p —)■ 0. Our 
dafa was more in favor of u{p) reaching a limif somewhere befween —0.6 and —0.8. In the light of 
[1] it was thus natural to revisit this earlier calculation and to derive the exact DSE for a link variable 
in Eandau gauge on the lattice. With this one could check if u{p) also on the lattice saturates the 
DSE in the infrared limit. Here we report on our first findings of a sfill on-going invesfigafion. 

2. Dyson-Schwinger equation of a lattice link variable in Landau gauge 

Eor fhe following if is advanfageous fo use a nofafion in which each link variable is assigned 
ifs defining sifes x and y = x±fL, i.e., we use fhe nofafion Uxy G SU{n) rather than the usual Uxp.- 
The Wilson gauge action for an SU{n) lattice gauge theory then reads 

1 ^ 

L with Pijki = i:v{UijUjkUMUu) (2.1) 

where E is a plaquette variable and g^ = 2n/^. Since we are interested in Eandau gauge we will 
consider links which minimize the (real-valued) Morse potential 

1 ^ 

nU] = -xY.^mj. ( 2 . 2 ) 

cj 

These minima fulfill fhe Eandau gauge condition (we use anfi-hermifian generafors U of SU{n)) 

0 = ff = Y.2IKtTv{fUj) ffa. (2.3) 

j 

Sfarfing wifh a fhermalized (non-gauged) configurafion U, gauge-fixing on fhe laffice is commonly 
performed by an iferafive procedure which consecutively gauge-lransforms U ^Ufj = giUijg^j until 
Eq. (2.3) is satisfied fo numerical precision. 

Eor the derivation of the DSE we now define an infinifesimal leff-variafion of such a 
gauge-fixed configurafion U = {Urs}' 

Sl^Urx = t%,„5rl5x^ - Un/5rmd,l + . (2.4) 
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It generates a left-variation of the link variable Uim followed by an infinitesimal gauge transfor¬ 
mation that returns the configuration to lattice Landau gauge (LLG). The anti-hermitian traceless 
nxn matrices {0,} in (2.4) ensure that the configuration remains in LLG. That is, the left-variation 
of the gauge condition is zero: 


0 = 6U—^HeTr 


t‘^t’’UU5u - 5i,„) + {diUij - Uijdj 


and the components of 0,- = solve the linear system, 

C,j 

where M is the Faddeev-Popov (FP) matrix 

k 


given in terms of the real symmetric matrix, 

^ab ^ ^ba rpj. . 


(2.5) 


(2.6) 


(2.7) 


The Faddeev-Popov matrix Mfj’ = My? is symmetric under the simultaneous exchange of color 
and site indices when Eq. (2.3) holds. Global gauge invariance of V implies that = 0, i.e., 

Mfj’ has rp- — I generic zero-modes. The remaining (A^ — l)(n^ — 1) eigenvalues are positive at 
each local minimum of V, and Eq. (2.6) implies that one may choose the solution 0“)^ of (2.5) 
orthogonal to the zero modes. 

If no gauge was fixed, the integration measure would be the Wilson measure djXy/ = 77 [L] . 

D\U] is the product of Haar measures dUij for each oriented link variable which is invariant under 
left- as well as right- group multiplication, d{gUij) = d{Uijg) = dUij for any g G SU{n), whereas 
Svr is invariant only under lattice gauge transformations Ufj = giUijg^^. 

To account for (minimal) Eandau gauge we introduce (in a Eaddeev-Popov-like manner) a 
density 

Pa([/) with , (2.8) 

j 


which in the limit a —)• oo has support only at the absolute minima U of Sic and whose integral over 
the gauge orbit jWidgi Pa{U^) = 1. In a first attempt we set Slg[U] = V[U], which for a — ^ oo 
gives support only for U, the global minimum of V on the gauge orbit of U. We will later see that 
for any of the standard lattice implementations of minimal Eandau gauge (these typically find only 
local minima of V) the identification =V[U] accounts only for the leading contribution in 

the DSE.^ Anyhow in the limit a —)• oo it gives for ‘2sa [67] 

llim [U] = ^detMp] (if Slg^V) (2.9) 

*The effective form of was unclear to us at the time of the conference, but it will be further discussed and 

specified in a forthcoming publication. 
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and the expectation value of a gauge-variant quantity 0 in minimal Landau gauge (MLG) becomes 


MLG = lim 7^ 

a^OO 


^ \d\U\ Pa(U} 
'W J 


( 2 . 10 ) 


with Zy/ = f D\U]pa,{U) . Note that in lattice perturbation theory one commonly sets Slg [U ] = 

ffff which gives |detM[17]| for 3!a[U] (see, e.g.,[4]). 

The Dyson-Schwinger equation (DSE) of a link variable in minimal Landau gauge we now 
obtain from a left-variation (see Eq. (2.4)) of the expectation value {Urs)MLG- expectation 
value must not change under the left-variation and so 

Applying to each term gives the DSE in an implicit form 

lim UlUrs + UrALiSw[U]+lnpa[U])) =0. (2.12) 

a-S’oo \ / MLG 


which after some algebra becomes 


Cf {TvUlm)MLG Wsm - 5„n5,/) = ^zY.^Kf,„TvfU, 


rs/MLG + IL 

ab 


na-,b a)/ba _ naM ot/baX 
^r,lm ^rs ^s;lm ^sr j 


MLG 


(2.13) 


Here Cf = {n^ — \)/{2n) is the quadratic Casimir invariant of the fundamental representation of 
SU{n), and the components of dj = Eq. (2.5). The (conserved) current is of 

the form 




Im 


Im * 


(2.14) 


where arises from varying the Wilson action, 


= \xRcY.Tr0U,,nUmjUjkUki , 

Sb j,k 


(2.15) 


and <I> from the variation of the induced measure Pa 

^UU] ■= lim 6Lln(pa[L]) = - 5u)Tr{t¥5ji-0t‘’5j,nyUi,n • (2.16) 

o;—>oo Z .Tt 

i j,bc 


Here we have used that Poo[U^] is stationary with respect to gauge transformations. This also 
implies that <I>"^ is transverse, = 0, in fact this is true at any minimum of V. Also is 

transverse in the sense that Z“^ = 0. 

The longitudinal part of the lattice DSE is algebraically satisfied by the last term of (2.13): 
Summing on the index "s" and using (2.3) and (2.6) we have at any minimum of E[f7] 

I " 5r„r)2l\z TrL/^ . (2.17) 

ab^s ab,s 
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3. Numerical verification 

To numerically verify our DSE we have Fourier-transformed Eq. (2.13) to momentum space. 
In this way we take advantage of its translation invariance to maximally reduce the statistical noise 
of the various Monte-Carlo expectation values. In momentum space the DSE reads^ 

=Z^v{p)+^^y{p)+Lf^y{p) (3.1) 

where y = ly) jin i?, the (momentum-independent) expectation value of the Morse potential in 
minimal Eandau gauge, while the terms on the rhs. are momentum-dependent. To verify that these 
dependences exactly cancel in the sum we have calculated (on several gauge-fixed ensembles) the 
term {ug = — 1) 

LfAp) = ^ E - OytnAf) <3-2) 

ah,xy ^ ' 

with 0 solving Eq. (2.5) and the transverse terms 

-TifVyy) 

^0 a,xy a 

and 

= s!r I ■ 

aycy 

Wxg is the sum of staples attached to the link and <I>"^ can be read off from Eq. (2.16) iden¬ 
tifying I = x and m = x + jX. For its evaluation we need different elements of [f/] which we 
estimate with the stochastic noise technique. It is remarkable, that a number of 8 to 32 Gaussian 
noise vectors for each U is sufficient to provide us with a good signal for <I>^^. 

In the continuum limit L^y becomes longitudinal. On the lattice this is not strictly the case, 
numerically however, the transverse contribution is negligible at the considered values of j3. 

In Fig. I we show the longitudinal (left panel) and transverse (right panel) terms of the DSE 
in momentum space. To this end, we have projected all terms with the respective longitudinal 
and transverse projectors, El = Pgpy/p^ and Pj = I—Pl where ap^ = 2sin(7rk^/L^). The figure 
hence shows the “form factors” of L^y, L^y and d>^v versus momentum a^p^. This is similar to 
what one typically does for the gluon propagator. 

In the left panel of Fig. 1 one clearly sees that the longitudinal part of the DSE is fulfilled for 
all momenta, as expected. Only L^y{p) gives a contribution, and the longitudinal part of L^y{p) — 
y 5^y = 0 (see the small panel on top). 

The transverse channel of the DSE has contributions from L^y{p), ^^y{p) and L^y. Their 
sum however does not equal yS^y. In fact, we see a clear deviation from zero for the difference 
'^S^y — (r^v + Lgv) which can be compensated if one rescales <I> by 1.3 (see top right 

panel of Fig.l). Without rescaling the deviation is proportional to the gluon propagator D^y{p), 
calculated for the same lattice parameters. This we verified for SU(2) and SU(3) for different j8 
and lattice sizes (see Fig.2). That is, the DSE would be fulfilled if if included anofher term cD^y. 

^Note, in what follows we change our notation for the links back to the usual, i.e., Uxy —> Uxfi- 


(3.3) 

(3.4) 
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5(7(3) : /3 = 6.0, 32* 5(7(3) : f = 6.0, 32* 




Figure 1: Longitudinal (left) and transverse (right) terms of the DSE as a function of the lattice momentum. 
The small panels on top show the validity of the (modihed) DSE for all momenta. 


5(7(3) : /3 = 6.0, 32* 
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5(7(2) : l3 = 2.3, 56'‘ 



Figure 2: Deviation of the transverse part of our initial ansatz for the DSE (full squares) in comparison to 
the corresponding gluon propagator times a constant (open diamonds). All in units of the lattice spacing and 
versus the lattice momentum squared. Left; SU(3) at j3 = 6.0,32^; Right; SU(2) j3 = 2.3,56^. 


At j8 = 6.0 the proportionality constant c is about 0.035. Our current ansatz for the DSE thus does 
not fully account for the Monte Carlo data at small momenta. 

We have not yet found a full explanation for this deviation. A subsequent analysis suggests that 
the reason for the deviation is our ansatz for lattice Landau gauge. Above we set Slg = V but our 
data at finite j3 effectively appears to favor local minima of Slg = V + £/^- Such an effective Morse 
potential would have the same minima as V but a different Hessian and for small e this Hessian 
would result in an additional term in the DSE proportional to the gluon propagator. At j3 = 0, 
however, the deviation is not longer proportional to the gluon propagator. At jS = 0, the rescaling 
<I> —)• 2<I> would approximately restore the DSE for the considered range of momenta (12"^). Eor 
more details, and hopefully an explanation, we have to refer to a forthcoming publication. 
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4. Kugo-Ojima 


Besides verifying the DSE we also want to check if the Kugo-Ojima correlator satisfies our 
DSE at p = 0. In the continuum the Kugo-Ojima correlator is given by the expectation value 
'Lb { (^vc)^(D^c)^) . On the lattice this corresponds to > 

where {c°Cy)cc = (^ 0^ is the Eaddeev-Popov matrix. We look at this correlator again in mo¬ 
mentum space where it reads 


xy abc ' 

Since we are interested in the limit u{p^ —)• 0), 
we have calculated u{p^) for the case of SU{2) on 
a 56^ lattice at j8 = 2.3. This has allowed us to 
reach relatively low momenta. This data is shown 
in Eig. 3, and one clearly sees the momentum de¬ 
pendence behaves as expected: u<x p^ for p^ —)■ 0. 
Nonetheless the limit u{p^ —)• 0) does not equal Y, 
shown as dashed line in Eig. 3. If this is a feature 
signaling the non-applicability of the KO criterion 
for minimal lattice Eandau gauge, or related to our 
ansatz for the DSE needs to be clarified yef. 

5. Summary 


y+v ^yv \ L ^x\i ^x~\-fi ^xfi J / * 

0.5 


0.4 

_ 0.3 
^ 0.2 


0.1 

0 

0 1 2 3 4 5 6 

2 2 
a p 

Figure 3: KO function vs. (lattice) momentum 
squared on a 56^ lattice; j3 = 2.3 (n = 2). The 
dashed line marks the value of Y. 



We have developed, for the first time, the Dyson-Schwinger equation for a lattice link variable 
in minimal Eandau gauge. The longitudinal channel of our DSE is algebraically satisfied, but for 
the transverse channel we see clear deviations whose origin is not fully understood, but will be 
further analyzed in a forthcoming publication. Once this is settled, our data for the KO function 
will also be revisited again. At present it signals the non-applicability of the KO criterion for lattice 
gauge theories in Eandau gauge, but this may be related to our present DSE. 
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